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Abstract 

We initiate a systematic study of the state space of non-extremal, stationary black hole 
bound states in four-dimensional N = 2 supergravity. Specifically, we show that an exponential 
multitude of classically stable "halo" bound states can be formed between large finite temper- 
ature D4-D0 black hole cores and much smaller, arbitrarily charged black holes at the same 
temperature. We map out in full the regions of existence for thermodynamically stable and 
metastable bound states in terms of the core's charges and temperature, as well as the region 
of stability of the core itself. Several features of these systems, such as a macroscopic config- 
urational entropy and exponential relaxation timescales, are similar to those of the extended 
family of glasses. We draw parallels between the two with a view toward understanding complex 
systems in fundamental physics. 
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1 Introduction 
1.1 Hot Halos 

Multicentered black hole and more general multihorizon solutions to Einstein's equations have been 
an active research area for almost a century [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 18, 19, 
20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37]. The first solutions discovered 
had string-like singularities connecting the holes. Others were time dependent or required highly 
fine-tuned exact cancelation of forces allowing arbitrary superpositions. Stationary, strictly stable, 
nonsingular solutions in four dimensions were only discovered more recently [11, 12, 14]. These 
solutions can have an arbitrary number of centers and represent genuine, molecule-like bound states, 
with equilibrium positions determined by balancing conditions between electrostatic, gravitational 
and scalar forces. 

Supersymmetry was instrumental in the discovery of these black molecules, but is in no way 
essential for their existence or stability. Indeed, since the black holes are effectively bound together 
by potentials with deep minima, deforming the supersymmetric solutions slightly away from ex- 
tremality, say by throwing some neutral particles into the black holes, will not destroy the bound 
state. Of course, if one keeps on adding mass to the system, the gravitational pull will eventually 
overpower all the other forces, causing the black molecule to collapse into a single black hole. Nev- 
ertheless, for sufficiently cold systems, one expects a plethora of stationary bound states very much 
like in the supersymmetric case. One of the goals of this paper is to make this argument more 
precise, and to map out part of the parameter space where such bound states must exist. 

We will do this in the setting of 4d N = 2 supergravity, by considering finite size black holes 
in the presence of a giant "galactic" black hole, all at the same temperature, taking the giant hole 
large enough so the finite size ones can be treated as probes in a fixed background — essentially a 
finite temperature version of the setup of [38]. This is rendered (artistically) in fig. 1.1. We argue 
that in this limit the probes can be taken to be BPS without significant loss of generality. Aided by 
the existence in the theories at hand of nontrivial families of non-extremal background solutions, 
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Figure 1.1: Left: A single probe is allowed to move on a sphere with radius depending on the relative 
values of its charges. Right: Multi-probe bound state configurations consisting of many probe halos 
orbiting a giant "galactic" black hole. To lowest order in the probe approximation the mutual 
interactions between the probes are neglected. The lowest energy bound states are stationary: No 
probes are actually moving, but there is angular momentum stored in the electromagnetic fields, 
pointing along the radial direction. Besides spin, the probes induce electric and magnetic dipole 
charges. In a holographic setup (e.g. embedded in AdS4) these would correspond to inhomogeneities 
in charge densities and magnetic fields in the dual field theory on R x S 2 . 

an exact computation of the probe interaction potential for any charge is accomplished. It is then a 
straightforward matter to determine the existence of local minima of this potential and to map out 
the subset of charges and temperatures where bound states exist and where they are energetically 
stable. 

We find that many of the interesting features that made their supersymmetric relatives rich 
and famous, such as wall crossing [51], field-induced spin [11], intrinsic Landau level degeneracies 
[13] and so on, persist at finite temperature. Close to non-BPS extremal limits, we observe that 
there always exist huge numbers of energetically stable probe bound states, that is to say bound 
states with free energy less than a single black hole with the same temperature and charge. Close 
to BPS extremal limits on the other hand, bound states tend to have slightly more energy than the 
black hole, so although classically stable, they can quantum tunnel into the horizon (though never 
to spatial infinity). But even on this side we find a region of parameter space where stable bound 
states exist upon turning on a temperature. 

We hope these results will facilitate searches for exact solutions of nonextremal black hole bound 
states, or for indirect criteria guaranteeing existence and stability. 
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1.2 Galactic Glasses 

Our primary motivation for this work, however, is to further the program of understanding complex 
systems in string theory in general and particularly, the holography of glass-like systems (other 
approaches to holographic glasses were proposed in [39, 40]). By complex or glass-like systems we 
roughly mean systems which have free energy landscapes with an exponentially large number of 
minima. Typically this leads to exponentially slow relaxation, breaking of ergodicity and persistent 
memory effects (see for example [41]). Examples in the real world include spin glasses [42], i.e. spins 
randomly distributed in a rigid (quenched) environment, effectively giving rise to a disordered mix 
of ferromagnetic and antiferromagnetic spin-spin interactions. Another and much larger class of 
examples are structural glasses, like glass itself, or electron glasses [43, 44], vortex glasses [45] or 
ice cream [46]. Here the disorder does not descend from randomly quenched auxiliary variables, 
but is generated spontaneously by the system itself, at least when cooled sufficiently fast. Much 
progress has been made in understanding the nature of glass and the glass transition over the past 
thirty years, yet it remains one of the deepest unsolved problems in solid state physics [47, 41]. 

Other systems sharing these general glassy features include biological systems such as the brain 
and other highly complex interacting networks. In string theory and cosmology, similarly complex 
systems arise in the description of the landscape, in weak coupling descriptions of black hole mi- 
crostates obtained from wrapped D-branes, and in the context of eternally inflating spacetimes. 
Some background material for this was given in [48, 49]. 

The basic observation suggesting black molecules have a holographic interpretation as thermo- 
dynamic states of glassy systems is quite simple. On the one hand, black holes can be thought of as 
thermodynamic states (in a precise, microscopic way in the context of holography). On the other 
hand, whenever black molecules exist, there exist exponentially many of them, including highly 
complex configurations with many centers. In different regimes, macroscopically different config- 
urations can be thermodynamically dominant [51, 71, 36]. Thus, a holographic dual description 
should have microscopic physics leading to a free energy landscape with exponentially many local 
minima and high complexity. Indeed, already at weak coupling, wrapped D-brane systems do have 
such energy landscapes [50, 51, 48]. 

Furthermore, as mentioned earlier, black molecules will eventually collapse when heated up or 
otherwise brought far from extremality. This process is classical and fast, and can be thought of as 
the equivalent of glass melting and entering the liquid phase. Going in the other direction however, 
even if thermodynamically favored, involves exponentially large relaxation times. Starting with a 
black molecular configuration that has a charge distribution over the centers that is macroscopically 
different from the minimal free energy configuration, the time scale to get to the equilibrium state 
may easily become exponentially long, for several reasons. First, the escape rate of a particle from 
a black hole is proportional to e~^ s , where AS is the amount of Bekenstein-Hawking entropy the 
emission costs the black hole [52, 54, 55]. Since charge needs to be moved from one black hole center 
to another to reach equilibrium, AS generically has a gap, due to charge quantization and the BPS 
bound. Since entropy scales quadratically under uniform rescaling of charge and mass, a generic 
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AS for the emission of a charged particle will scale linearly with the size of the black hole, leading 
to typical relaxation time scales that are exponential in the system size. Second, even if there 
happen to be some relatively fast decay channels away from a particular configuration, relaxation 
is likely to get slowed down tremendously by the existence of an exponentially large landscape of 
black molecule configurations, much like what happens in actual glasses. Indeed there is in general 
no reason why the fastest decay channels out of a particular configuration should drive the system 
straight to the configuration that globally minimizes the free energy. Instead the system may get 
stuck in deep local minima of the free energy, requiring again exponentially suppressed transitions 
to get unstuck. 

Finally, it seems plausible that these systems exhibit memory effects similar to those observed 
in glasses. When a glassy material is perturbed by a constant force for some time t p and then 
allowed to relax back for a time t, one frequently observes that it approaches its original state as 
5(t) ~ log(l + tp/t) over many decades of time [43, 44]. A striking feature of this is its completely 
universal form and the absence of any intrinsic time scale, replaced instead by a persistent "memory" 
of t p . Such behavior is explained robustly [43, 44] when one assumes a high dimensional state 
space with linear stochastic dynamics given by exponentially small transition matrix elements 
Tij ~ e y , and roughly uniformly (say polynomially) distributed values of Aij. The density 
of transition channels with relaxation time T = e A is then proportional to ^ at large T, i.e. 
approximately scale invariant. This naturally leads to an approximately scale invariant spectrum 
of the transition rate matrix, producing the above logarithmic time evolution. Since the systems 
of interest to us exhibit a similarly broad quasi continuum of exponentially suppressed transition 
rates, it is reasonable to expect similar memory effects to arise. 

Thus it is conceivable that these systems are part of the large family of glassy systems. On the 
other hand, it is far from obvious if and to what extent they behave like conventional glasses, what 
their properties are and what their holographic interpretation is, or what the potential lessons are 
we can draw from their study regarding more general complex systems in string theory. In this 
paper we will take the necessary first steps to address these questions. 

1.3 Layout of the paper 

The main focus of our work consists of exploring the nature of finite temperature black hole bound 
states in N = 2 supergravity in four dimensions. To be more or less self-contained we begin by 
reviewing in section 2 the relevant Lagrangians and notations (more background can be found in 
[56, 57]). In section 3 we present a simple consistent truncation scheme valid for any J\f = 2 string 
compactification, within which we derive and describe a family of exact non-extremal black hole 
solutions (previously found in [58]). These will serve as our galactic black hole cores, dressed by 
halos of much smaller probe black holes. The probe potentials are introduced in section 4 and the 
existence of nonextremal bound states is established. A systematic exploration of the parameter 
regime in which metastable and stable bound states exist is given in section 5; the results are shown 
in fig. 5.1. We discuss the thermodynamic properties of the system and find a phase structure 
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confirming several of the glass-like features outlined above. The diagram is suggestive of quantum 
critical points attained by dialing the asymptotic moduli to the black hole attractor fixed point. 
Across this point, the ratio of probe-induced spin to probe induced D6 magnetic dipole moment 
changes sign. In the appendix we estimate the configurational entropy of multi-probe BPS galaxies 
and find it grows linearly with the system's size (charge) , and that it scales to zero with a nontrivial 
exponent near the critical points. 

Note added: In interesting recent work [89, 90], analogous but complementary finite temperature 
bound states in five dimensions were independently explored, with qualitatively similar results (see 
appendix D of [90] for a detailed comparison). 

2 Setup and notation 

Four dimensional N = 2 supergravity coupled to massless vector and hypermultiplets has a bosonic 
action of the general form [59, 60] 

S 4D = ^ I d 4 x V=g (\R - Gab d^ A d»z B - h XY d^q x d^q Y ) 
+ lL J d^^g^ATjjF^F^ -ReAfjjF^F^) , 

where the z A {A = 1, . . . , n) are the vectormultiplet scalars, F l = d^A u — dyA 1 ^ [I = 0, 1, . . . , n) 
are the vector field strengths, F^ v = \e^ upa F pa ', and the q x are the hypermultiplet scalars. We 
put Gn = 1- As G A g and A/}j only depend on the z A , and hxY depends only on the q x , 
the hypermultiplets decouple from the vector multiplets and we will not need to consider them 
further. The vectors A are sourced by electric charges Qj and magnetic charges P 1 . l The space 
of charge vectors V = (P^Qj) carries a canonical, duality invariant, symplectic product, which in 
the standard symplectic basis can be expressed as 

(r,f) =P I Q I -Q I P I . (2.1) 

The metric G A g is special Kahler, i.e. it is derived from a prepotential F{X): 

G A B = d zA d- S lC, /C = -logi(«,n), n = (X 1 ,d x iF) , X A = X°z A . (2.2) 

The variable X° drops out of all observable quantities; we gauge fix X° = 1. The prepotential F(X) 
is a locally defined holomorphic function, homogeneous of degree 2 in the X 1 . It also determines 

1 Magnetic charges have an upper index, but we will further on often use lower indices for both electric and 
magnetic charges, to make expressions involving powers of magnetic charges less clumsy. 
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the electromagnetic couplings Nij: 

tfu-Fu + 2i X M {lmFMN)X N > *>J-^W- (2-3) 

In type IIA Calabi-Yau compactifications the coordinates z A = B A +iJ A are identified complex- 
ified Kahler moduli and, ignoring string worldsheet instanton corrections, the prepotential takes 
the form 

F(X) = -L Dabc X a X b X c (2.4) 

where the Dabc are triple 4-cycle intersection numbers. The charges (Pq, Pa, Qo, Qa) are identified 
with wrapped (D6,D4,D0,D2) charges. For a Calabi-Yau with a single Kahler modulus z = z 1 (i.e. 
n = 1), the cubic prepotential (2.4) becomes 

(X 1 ) 3 

FiX) = D ~6X^- (2 ' 5) 

The examples we consider in this paper will all be effectively reducible to this case. The special 
Kahler metric (2.2) is then just the Poincare metric on the upper half plane: 

3 

Gz ~ z = I(hn~zf ' (2 - 6) 
and the electromagnetic coupling matrix (2.3) becomes, with z = x + iy, 

( x 3 i iyx 2 , iy^_ _x^_ _ iyx\ 
3 + 2 2 + 6 2 2 (2.7) 
_x± _ ryx , iy I v > 

2 2 x ' 2 / 

N = 2 supersymmetry implies that the mass of any state of charge T in a vacuum with asymptotic 
moduli zq is bounded below by the absolute value of the central charge Z(F, z$), defined by 

Z(T,z) = -e K/2 {T,Q). (2.8) 

States saturating this bound are supersymmetric and called BPS. In the n = 1 case (2.5), for a 
charge V = (Pq, Pi, Qo, Qi), (2.8) becomes 



V3_ 
2 y/DfJmzY 



Z(T, z) = - y^== (f Poz 3 - f Piz 2 + Q x z + Qo) • (2.9) 



To write down the action of a point particle in a general background, it is convenient to introduce 
the dual magnetic field strengths 

Gi = Ira Nij F j - Re M u F J . (2.10) 
The electromagnetic equation of motion dGj = implies the existence of dual magnetic gauge 
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potentials Bj such that Gj^y = d^Bj^ — d u Bj u . Collecting the electric and magnetic gauge field 
strengths and potentials into duality covariant vectors F = (F 1 , Gj) and A = (A , Bj), the action 
for a point particle of mass m and charge 7 = (p 1 , qj) is [11, 61] 

S 1 = - [ mds- - f (7,A M )dx". (2.11) 



2 „ 

The mass depends on the scalars z. In particular when the particle is BPS, we have 

m = \Z(r/,z)\. (2.12) 

3 Non-extremal black hole background 

We will now construct a class of exact, spherically symmetric, nonextremal single centered black 
hole solutions for any prepotential of the form (2.4). These are essentially the solutions found in 
[58, 62]. In the nonsupersymmetric extremal limit they belong to the class studied in [63, 64, 65, 66]. 

3.1 Equations of motion 

The black hole metric is of the general form 

ds * = _ e 2U(r) dt 2 + e -2U(r) ( _A^ dT 1 + , (3.1) 

\sinh cr sinh cr / 

where r is an (inverse) radial coordinate with r = corresponding to spatial infinity and r = 00 to 
the horizon. The parameter c is a positive constant parametrizing the deviation from extremality. 
The scalars depend on r only, and the electromagnetic field is given by 2 

F 1 = Piuj - (lmJ\f) IJ (Qj + KqMjkPk) *w , u = sin 9 dO A d<p , *w = e 2U dt A dr (3.2) 

Note that this automatically satisfies the Bianchi identity dF 1 = 0, because the Mij and U depend 
on r. Moreover, for this particular form of F 1 , we have Gj = Qjoj + (••■) *w, hence the equations 
of motion dGi = are also automatically satisfied. 

The scalar and metric equations of motion can be obtained from an effective particle action [67] 

SW = dr (U 2 + G AB z A ^ - V cS (U, z, z)) , (3.3) 

with effective potential 

V cS = -c 2 - e 2U {\Z\ 2 + AG AB d A \Z\ B B \Z\) , (3.4) 



2 In form notation, F = ^F^dx^ A dx v , A = A^dx 1 * 
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supplemented with the constraint that the total particle energy must vanish: 

U 2 + G AB z A z- B + V eS = . (3.5) 

3.2 Consistent truncations 

Solving this system in general appears intractable, but special classes of solutions can nevertheless 
be found. First, the general problem with an arbitrary number n of vector multiplets can be 
consistently reduced to an effective single vector multiplet problem by the truncation [51], for any 
choice of constant K A (inside the Kahler cone): 

z A = K A z\ (F°,F a ,G ,G a ) = (f ,K a F 1 ,K 3 G ,(K 2 ) a G 1 ) , (3.6) 

where (K 2 ) A = D A bcK b K c and K 3 = D AB cK A K B K c . It is easily checked that the equations 
of motion then consistently reduce to the n = 1 case (2.5), with D = K 3 . To remain consistent we 
must also choose the black hole charge to be of the form 

T = (p°, K A P\K 3 Q , (K 2 )aQi) , (3.7) 

which sources the reduced fields as a charge (P ,Qi) in the effective n = 1 theory. We will 
henceforth normalize K A such that D = K 3 = 1, and drop the hats on the reduced quantities. 

The effective metric on the scalar space parametrized by z 1 = z = x + iy is given by (2.6) and 
the effective potential is V c g(U,x,y) = c 2 — e 2U V(x,y) with 

V = \ (Qo + & x - \ P ^ 2 + S P ox 3 ) 2 + -{Qi- Pix + \Pox 2 f + | (Pi - P x) 2 + ^P 2 . (3.8) 

The resulting equations of motion are still hard to solve in the generic case, but when x = 0, the 
coupling matrix (2.7) becomes pure imaginary and diagonal, and the system simplifies considerably. 
This motivates a search for solutions with constant x(r) = 0. In this case the fields strengths (3.2) 
are of the simple form 

u ^u y a > i *°6i y i (3.9) 

G = Qqoj + Po\ *oj , G\ = Qi oj + Pi | *u) . 

Consistency of the ansatz x(t) = requires d x V\ x= o = 0, which leads to the conditions QoQi = 
QiPi = P Pi = 0, leaving the possibility to have D4-D0 (r = (0, Pi, Q , 0)), D6-D2 (r = 
(Po, 0, 0, Qx)), or D6-D0 (r = (Po, 0, Qo, 0)) background charges. In the following we specialize 
to these cases. 
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3.3 Solving the equations of motion 

Putting y = e^, the remaining equations of motion for U and <fi derived from (3.3) take the form 



U = e 2U ( Vl e a(t> + v 2 e b +), (3.10) 
4> = \e 2V (av^ + bv2e H ) , (3.11) 

together with the constraint (3.5) (which just fixes the value of c). For the D4-D0 system, we 
have (v 1 ,v 2 ;a,b) = (zQl, ^ ; -3, l) , for D6-D2 (vi,v 2 ;a,b) = (q%,$;-1,£) and for D6-D0 

(vi,V2',a,b) = ^3Qq, j^; —3,3^. This system is of Toda form [68]. Following the method of [62], 
we set a = 2U + a<f>, (3 = 2U + b<p. The system of equations of motion for U and <j) then becomes 

q = a e a + 7 e /3 , 
P = 5 e a + f3 oe P. 

where ao = § (3 + a 2 ) v\, /3q = | (3 + b 2 ) v 2 , 7o = § (3 + ab) v%, and Jo = § (3 + ab) v 2 . These two 
equations decouple if ab = —3, which happens to be the case for the D4-D0 and D6-D2 systems. 
In these cases the equations of motion integrate to 



a (r) = log 
Pir) = log 



2c? 



«o sinh 2 (c]T + c 2 ) 



2c\ 



i sinh (c 3 r + c 4 ) 



(3.12) 



where ci, 02,03,04 are positive integration constants and (ao,A)) = ^24<3g, for the D4-D0 

while (q ,/3 ) = ^) for the D6-D2. 

Specializing to the D4-D0 case, this implies for the original fields 



e -AU = 2\QoPl I ginh , + ginh 3 (cgr + , (313) 

3ClCg 



y 2 = e 2<f> = 6|Qp| c 3 sinh (cir + c 2 ) 
|Pl| ci sinh (c3T + C4) 

The constraint (3.5) fixes c 2 = (c 2 + 3c 2 )/4. Regularity of <f> at the horizon r = 00 requires ci = C3. 
The asymptotic boundary conditions £7(r = 0) = 0, y(r = 0) = yo further imply 

sinhc2 = 7T75^ 1 ' sinhc 4 = , ^1/0 • ( 3 - 15 ) 
2V3 |Qo| \Pl\Vq 
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3.4 The D4-D0 solution 

Putting everything together, denoting 3 



Hq = sinhfcr + c 2 ) , H\ = - — — sinhfcr + C4) , (3.16) 

c c 



we get for the metric warp factor and the scalar 



e~ 2U 



\h*HI (3.17) 



</ - A/^, (3-18) 



and for the gauge potentials A = (A 1 , Bj), obtained by integrating the field strengths (3.9): 

A ° = Ao {\j° 2 + If " c ) dt ' A 1 = P 1 (l-cos^)#, 

/ / \ (3.19) 

S = Qo(l-cos5)# ; B i = -wi (Jc 2 + hJ 



We have chosen a gauge here in which the electric potentials vanish at the horizon and the Dirac 
monopole potentials are regular on the northern sphere. 

Notice that the modulus y at the horizon r = 00 is fixed at an attractor point independent 
of yo and c: 

9 *=/if- (3 - 20) 

3.5 Mass, entropy, temperature and specific heat 

The ADM mass of the black hole can be read off from the asymptotic behavior of the metric: 

M = |(cothc 2 + 3cothc 4 ) (3.21) 




c2 + ^j_yo _ (3 22) 

In the extremal limit c = 0, this becomes M = -^l-Pilyo^ + ^-^i^- ^l 1611 ^ an( ^ Qo have the same 
sign, this equals the absolute value of the central charge and the extremal limit is supersymmetric. 
When Pi and Qo have opposite sign the mass is strictly larger than the absolute value of the central 
charge and the extremal limit is nonsupersymmetric. 

The Bekenstein-Hawking entropy can be read off from the near horizon (r — > 00) behavior of 



3 This notation is motivated by the fact that in the c — > extremal limit, Ho and Hi become the flat space DO 
resp. D4 harmonic functions ubiquitous in the description of the well-known extremal solutions. 
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Figure 3.1: (a): Rescaled temperature T = ^/yo~Pi T as a function of Qq and c, multiplied by 10 3 . 
Red is warm, blue is cold. The temperature reaches its maximum at fixed Qq on the thick curve, 
(b): Specific heat for Qq = 0, diverging at c cr i t = l/2y/6 ~ 0.2, where the temperature reaches a 
maximum T max = V3/8tt sa 69 x 10-3. 



the metric. Introducing r = e CT , this is: 



S vr 

where S = ^4hor/4 is the Bekenstein-Hawking entropy: 



r 2 dt 2 + —dr 2 + -dn 2 



TT 



2 > 



(3.23) 



'2\QoP?\ 



,(c 2 +3c 4 )/2 



vr c + 




12 Qg 



1/2 



C+ \/c^ + 



\ 3/2 



(3.24) 
(3.25) 



The temperature can be read off from the near-horizon metric by Wick rotating the time 
coordinate and fixing its periodicity f3 = 1/T by requiring regularity at the origin r = 0. This 
yields 

T =2S- (3 - 26) 



Whereas the mass and entropy of the single centered D4-D0 black hole are monotonic function 
of c, this is not so for the temperature, as is clear from figure 3.1 (the rescaled tilde variables will be 
defined in (3.27) and below): The temperature starts at zero in the extremal limit c = 0, acquires a 
maximal value at some intermediate c, and goes to zero again in the Schwarzchild black hole limit 
c — > oo. The rescaled specific heat O^f at fixed Pi,Qo,yo diverges at this point [69] and changes 
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sign from positive to negative. 



3.6 Scaling symmetries 



It turns out to be very useful to keep in mind the scaling symmetries X 
acting on the various quantities defined so far with exponents 



\?\?X, Xi e 





Pi Qo c yo t 


Ho Hi 


M S T 


m 

n 2 


1 110-1 
1 3 | 1 -| 




3 1 
2 2 


12-1 

3 q 3 
2 2 



The first symmetry descends from a general symmetry of Einstein gravity coupled to vectors, the 
second from a general symmetry valid for cubic prepotentials [51]. The scalings of the derived 
quantities Ho, Hi, M, S and T follow from the scalings of the charges and of c, yo an d t. A conse- 
quence of these symmetries is that physical quantities will depend only on invariant combinations 
of the parameters, up to an overall factor determined by the scaling properties of the quantity 
under consideration. We choose our independent invariant parameters to be 



Qo 



Qo 



y 2 Pi VvoPi ' 

A quantity X with scaling exponents {n\,n2) will then have the functional dependence 
X{Pi,yo,Qo,c) = P^y^- ni X{Qo,c) , X(Q ,c) = X (1,1, Q ,c) . 



(3.27) 



(3.28) 



For example the ADM mass (3.22) can be written as M = Pi^/y^M with M = c 2 + 12 + 
|yc 2 + the entropy as S = PfyoS and the temperature as T = T/Pi^/y~o~. 



4 Probe bound states 

In this section we explicitly demonstrate the existence of bound states of probe particles in the 
nonextremal D4-D0 black hole backgrounds described in the previous section. 

4.1 BPS probes 

We will in this paper primarily consider probe particles of charge 7 = (po-,Pii <7o> Ql) that are 
themselves BPS. These "particles" could themselves be large black holes, but they must be much 
smaller than the background black hole so backreaction can be neglected. In thermal equilibrium 
the probe black hole will acquire the temperature of the background black hole, so it will not quite 
be BPS. However, the background temperature is parametrically suppressed as T ~ 1/M in the 
limit of large background black hole mass M, as can be seen explicitly from the scaling table in 
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section 3.6. Hence for any fixed probe size, the thermal contribution to the probe energy, which is 
proportional to T, will vanish in the large M limit, and thus the BPS approximation is justified. 4 
The static potential for a BPS particle can be obtained from (2.11) and the solutions found 
in the previous section. It consists of two parts, a gravitational part V g = e u \Z('j,y)\ and an 
electromagnetic part V cm = ^(7, Aq). Explicitly V p = V g + V em with 




and 




^ = -7^h ^ + ^ 2 -c -7^ A h 2 + ^k-c\ , (4.2) 



To avoid complications with marginal stability decays of the probe as it moves around in the 
nontrivial background, we will only consider probes that are themselves single centered black holes 
or particles. The BPS probe entropy is given by S p = irVT>, where the so-called discriminant T> 
must be positive for a solution to exist; in the case at hand this is [70]: 

V = \ ~ P0Q0 ~ 2 PoPiQoqi + ^phl - ^P°1i - °- ( 4 - 3 ) 
Finding an explicit parametrization of this subset of charges seems hard, but is actually made easy 

2 

by using the invariance of V under shifts (pojPii <Zo ; Qi) = 7 — > 7a = (po,Pi +Po a , Qo — 9l« ~ Pi% — 

3 2 

P073") Qi +Pifl+Po^") with a£l, This invariance follows from the axionic shift symmetry z — > z — a 
of the supergravity theory under consideration, which leaves in particular the black hole entropy 
invariant. 5 

Thus, in the parametrization 

k 2 k 3 
Pi/Po = k, qi/po = -b + — , qo/po = n + bk - — , (4.4) 

2 D 

we get simply 

V = pUlb 3 -n 2 ), (4.5) 

1 /3 

and we may explicitly parametrize the solutions to the constraint (4.3) as b = ^§("< 2 + , 
V > 0. 

In type II A compactifications, k may be thought of as the U(l) flux on the wrapped D6, which 
carries no entropy, while b and n are the "entropic" contributions to the charges [51]. 



4 More precisely, the probe thermal energy is of order Et ~ TS P , with S p the probe entropy. Scaling up the 
background charges uniformly by a factor of Ai while keeping the probe fixed scales Et oc T oc f/Ai — > 0, whereas 
the probe potential remains invariant (as we will confirm below). Assuming this potential is not exactly fiat, we 
can therefore neglect the thermal energy. Strictly speaking, this argument only tells us we can assume the probe to 
be extremal, but not necessarily BPS. However non-BPS extremal black holes are expected to be unstable to decay 
into lighter particles, as will be confirmed explicitly in section 5, and this on a time scale exponentially smaller than 
possible instabilities of the background black hole. This justifies considering primarily BPS probes. 

Explicitly: 0( 7 ) =min,|Z(7,z)| 2 = min 2 \Z(j, z - a)\ 2 = min, \Z(j a , z)\ 2 =£>( 7a ). 
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4.2 Scalings and validity of probe approximation 

Besides the scaling symmetries described in section 3.6, we have an additional symmetry uniformly 
scaling only the probe charge 7, present because we work to linear order in 7. All in all we get the 
following scalings X -»• X^X^X^X: 





Pi 


Qo 


c 


yo 


T 


Po 


Pi 


Ql 


qo 


k 


b 


n 


m p 


Sp 


v P 


m 


1 


1 


1 





— 1 
































n-i 


1 


3 


3 
2 


1 


3 
2 





1 


2 


3 


1 


2 


3 


3 
2 


3 




"3 

















1 


1 


1 


1 











1 


2 


1 



The action on the background is the same as before, and we used the third scaling to set to zero the 
action of the Ai-scaling symmetry on the probe charge. In addition to these continuous symmetries, 
there is a Z2 symmetry inverting the signs of all charges, and another Z2 inverting only the signs 
of D2 and D6-charges. We point out here that the probe potential does not scale with the size of 
the black hole. 

Analogous to (3.27) and (3.28), we can scale out powers po in addition to Pi, yo according to 
the following scaling dimensions: 

X = P^y^pfX. (4.6) 

So for example k = y^k, b = y^b, n = yjjn, and S p = y^p^Sp with S p = tt^J |& 3 — n 2 . 

Some care has to be taken not to forget the regime of validity of the probe approximation. 
The ratio probe mass m p over background black hole mass M has scaling weights (—1,0,1), so 

= yo£o _ F rom thi s we see that for typical tilde variables of order 1 (and po 7^ 0) , we must 
keep yopo <C Pi to guarantee m p /M -C 1. Indeed, when yo ~ Pi, a single pure D6 becomes as 
massive as P\ D4-branes, i.e. the background black hole, thus spoiling the probe approximation. 
In particular, since yo — > 00 is the M-theory decoupling limit [71], this means that black hole 
bound states in AdS3 x S 2 (of which the exact supersymmetric versions were constructed in [71]) 
are not reliably captured by the 4d probe analysis of this paper, in particular not for non-BPS 
configurations. On the other hand, for any fixed value of yo and the probe charges, we can always 
send P\ — >• co to make the probe approximation arbitrarily accurate. 

4.3 Bound states 

The probe will form a stationary "molecular" bound state with the black hole whenever the potential 
has a nontrivial local minimum. In the supersymmetric case, the discovery of such probe bound 
states led the way to the construction of general nonlinear black hole bound state solutions in M = 2 
supergravity [11]. In particular their existence made it clear that such bound states had to exist, and 
quite remarkably, the simple explicit formula for the equilibrium radius obtained from the probe 
analysis is formally exactly the same as the corresponding formula obtained from supergravity. 
This was later explained as being a consequence of the constraints imposed by supersymmetry [13]. 
There is no reason to expect a similar exact match in the nonextremal case, but a probe analysis 
will still provide reliable information about the existence of bound states in suitable regimes. 
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Figure 4.1: Examples of probe potentials. Local minima give rise to bound states. The radial u 
coordinate is denned by = g^og^ ■ The horizon is at u = and spatial infinity at u = 1. In 
general, the potential always goes up from the horizon (for c > 0) and has at most one bump and 
at most one interior local minimum. On the left we show probe potentials at (Qo, c, k, 9T>, h) equal 
to (.01, .01,. 5, 0,0) for (1), (.4, .01, .65, 0,0) for (2) and (-.4, .01, 2, .01, .15) for (3). On the right we 
show the same but at c = .08. At this higher temperature, the minima of (1) and (2) have become 
positive and more shallow, while (3) has lost its minimum altogether. Increasing c even more wipes 
out all local minima. 

To reproduce first the supersymmetric result, we consider the case of supersymmetric back- 
ground, c — > 0, Pi, Qo > 0. The probe potential is then of the form V p = -\Aem + ^ 2 + Km- 
This makes the BPS bound V p > manifest. If a BPS-saturating supersymmetric minimum 
V p = exists, it is reached at the radius r eq = l/r eq for which A(r eq ) = and V era (T cq ) < 0. 

3/2 
Vo 



The first condition is giPi(r cq ) — PqHq{t c ^) = 0, or explicitly, using lim^o Hq = Qqt + and 
lim c ^o#i = Pit + 

Vo 

PoQo - 91 Pi fA ^ 

This reproduces the standard BPS equilibrium separation formula for two-centered bound states 
of this kind [11]. 

In the nonextremal case there is no such simple expression for r eq , but by continuity there will 
obviously still exist bound states for suitable values of the charges and the nonextremality parameter 
c. Some examples of probe potentials with local minima are shown in fig. 4.1. As suggested by the 
figure, increasing the nonextremality parameter c typically tends to push up the local minimum, 
until it eventually disappears altogether and the probe rolls into the black hole. This is to be 
expected, since going away from extremality means adding more mass. Thus the gravitational 
pull becomes increasingly more important, eventually overpowering all other forces. In some cases 
however, in particular for small positive values of Qo, going away from c = initially decreases the 
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value of the potential at the local minimum. An example is potential (1) in the figure (as opposed 
to (2)): At c = this is a supersymmetric (i.e. V = 0) local minimum, whereas for small but 
nonzero c it is negative. For larger c is goes positive again. Thus, interestingly, slightly heating a 
supersymmetric black hole with small Qo will make it unstable to emission of such charges. We 
view this as an interesting interplay between supersymmetric and thermal physics. For non-BPS 
extremal black holes (Qo < 0), negative energy probe bound states exist for sufficiently small c for 
all values of Qo. An example is potential (3) in the figure. 

In the following section we discuss existence and stability in more detail. 



5 Existence, stability and phases 




Qo Qo 



Figure 5.1: Existence and stability regions of hot black molecules in the (Qo, c)-plane. The figure 
on the right zooms in on a smaller region but is otherwise the same as the figure on the left. No 
bound states exist in the grey region. In the yellow region bound states exists, but they are all 
positive energy (curve (2) in fig. 4.1). In the green region, negative energy bound states exist 
(curve (1) in fig. 4.1a). The orange overlay is the region where the black hole core itself is unstable 
for emission of particles to infinity. The blue overlay in the lower left corner is the region within 
the Qo < — g range where pure fluxed D6 probe bound states exist (when Qo > —\ these always 
exist and moreover they always produce the lowest energy bound states). The dotted line is the 
maximal temperature line also shown in fig. 3.1. The grey regions touch the zero temperature axis 
at Qo = igj or equivalently when yo coincides with the black hole fixed point y*. 
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5.1 Supersymmetric bound states 

In the supersymmetric case Qq > 0, c = 0, using the parametrization (4.4), a straightforward 
analysis shows that a BPS bound state exists if and only if 

Qo < 91 < \ or \ < qi < Qq , (5.1) 

D O 



where qi = \ — b, and 



sign k = sign po , V = -b — h > . (5.2) 



The equilibrium distance (4.7) expressed in rescaled variables is f eq = • ^he boundary 

P = corresponds to vanishing probe entropy, the boundary q\ = Qo to a vanishing bound state 
radius and hence absorption of the probe by the background black hole, and finally the boundary 
q~i = | corresponds to an infinite bound state radius and hence to decay at marginal stability. This 
is brought in a more conventional form by returning to the non-scaled variables, which turns the 
existence condition (5.1) into < < ^2 or 6<2o -> 621 ~> y2 The absorption wall is then 
clearly seen to correspond to a vanishing probe-background symplectic product (7, T) = 0, while 
the marginal stability wall is at yo = \J^~- 

Notice there exist bound states for all values of Qq > except 1/6. The number of possible 
bound states will not be constant however. In particular when Qo ~ ^ 1/6 (or equivalently yo ~~ * U+i 
where y+ = is the attractor fixed point (3.20) of the background black hole), the allowed region 
in the probe charge space shrinks to zero. In appendix A we compute the number of probe bound 
states, allowing multiple probes with different charges and taking into account the lowest Landau 
level degeneracies due to the magnetic interaction between the background black hole and the probe 
charges (but ignoring mutual magnetic interactions between the probes themselves). We do not 
count the internal microstates of the black holes. The logarithm of the number of configurations 
defined in this way is thus the analog of the notion of configurational entropy in the theory of 
glasses [72, 73, 74]. The final result for the number M(e) of such configurations with total probe 
mass over black hole mass less than e is given by equation (A. 18): 

iogAA( e )~6 5 / 6 | yo 2 -^r /3 4^' ( 5 - 3 ) 

with = ^p^-, as in (3.20). Thus we see that the number of allowed configurations indeed goes to 
zero when the critical point yo = U* is approached, with a nontrivial scaling exponent 1/3. 

It is interesting that even this restricted counting already gives an exponential growth of the 
number of configurations N in the P\ — > 00 thermodynamic limit. This is typical for glasses 
[72, 73, 74]. The growth is not as fast as the black hole entropy itself (it is at most P^ 2 , if we 
allow going to the boundary of the probe regime, cf. eq. (A.6)), but the exponentially large number 
of configurations should nevertheless have important consequences for the thermodynamics of this 
system. 
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5.2 Hot black molecules 

For nonsupersymmetric black holes the analysis becomes more complicated, requiring some numer- 
ical assistance to scan the space of possible probe bound states for given (Qo,c). The results of 
this work are summarized in fig. 5.1. We identify four different regions: 

1. In the grey region filling the high temperature region, no molecular bound states of any kind 
exist, as gravity overpowers all other forces. 

2. In the yellow regions right below it, bound states exist for some probe charges, but all of 
them have positive energy at their minimum, so they are metastable (recall we put the zero of 
the probe potential at the horizon). An example is potential (2) in fig. 4.1. When approaching 
the grey-yellow boundary, the minima become higher, are pushed to large radii and become very 
shallow, while the number of probe configurations goes to zero. This should give a scaling law 
analogous to (5.3) but we did not try to extract the scaling exponent. The grey region touches 
the T = axis at the quantum critical points Qq = ±g « 0.1667. For Qq > — g, we numerically 
observed with high accuracy that the probe particles forming the lowest energy bound states are 
always zero entropy b = n = 0, k ^ particles. In IIA language these are pure D6-branes with 
U(l) flux, which uplift in M-theory to smooth "bubbling" geometries [75, 76, 77, 78, 79, 80, 81]. In 
particular the bound state surviving the longest when c is increased is of this type. For Qq < — | 
this is no longer the case, and in fact there are no bound states of this type for — | < Qq < — g. 
The blue overlay shows where they reappear in the region Qq < — |. 

3. In the green regions enclosed in the yellow, negative energy bound states exist. Such bound 
states are energetically stable against tunneling of the probe into the black hole or out to infinity. 
The green line along the positive Qq axis represents the BPS bound states discussed in section 
5.1, which have zero energy. The negative Qq axis represents extremal nonsupersymmetric bound 
states. There is no BPS bound forbidding negative energy states, and by the rule that everything 
that is not forbidden is allowed, we find indeed that a large subset has deep negative energy minima. 
Interestingly, due to the transient dipping effect described at the end of 4.3, there is a small but 
finite region on the Qq > (i.e. BPS) side at finite temperature where negative energy bound states 
exist. It extends to Qq = ^ « 0.0185 (see end of next paragraph). 

4. The orange overlay is the region where the background black hole itself is unstable to emission 
of particles to spatial infinity. We take such emissions to be possible whenever there exists some 
probe charge such that the probe potential becomes negative at spatial infinity. Notice that the 
red region includes the green region. Hence whenever a bound state exists that is stable against 
tunneling of the probe out of its minimum, the background black hole will be unstable to emission 
of particles. We note, however, that whenever the probe potential exhibits a minimum, it is always 
found to be lower than the value of the potential at infinity. The destabilizing probe type kicking 
in first (at the highest c) is again a pure fluxed D6 brane, b = n = 0, over the entire range of Qq 
we scanned. This allows computing the boundary of the red region analytically as the value of c 
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Figure 5.2: Same as fig. 5.1, but for the system confined in a box, which can be viewed as a rough 
model for AdS. The box is implemented as a cutoff at radius sinh 6 (g f ) = 1, i-e. u = 1/2 (with u 
defined in fig. 4.1). Bound states with equilibrium positions at u > 1/2 are discarded, and decay 
is defined as emission to u = 1/2. 

for which the asymptotic value of V p has a double zero viewed as a function of k: 



The critical line reaches zero at A = g, i.e. Qo = 34) which is numerically seen to coincide with the 
edge of the green region. It asymptotes for Qo — > —00 to ^ w 0.433. 

5.3 Bound states in a box 

The black holes under consideration live in asymptotically flat space. Non-BPS black holes are 
unstable due to Hawking radiation and possibly other emission instabilities, and hot flat space 
is unstable due to nucleation of black holes [82, 83]. As a result, it is hard to make sense of 
this infinitely large system as a statistical mechanical model. To make it better defined, we can 
put the system in a finite box, either by imposing a cutoff by hand at some finite radius or by 
embedding the system in AdS4, along the lines of [84, 85, 86, 87], allowing the black hole to achieve 
thermal equilibrium with its environment in the box. In simple setups that do not give rise to 
black molecules, one sees that depending on the size of the box and the temperature, the final 
equilibrium state can either be a big black hole , or an ordinary thermal gas. 6 In the case at hand, 
we may expect much more complicated, glass-like behavior at sufficiently low temperatures, due to 
the exponentially large number of complex stationary configurations that exist. 

6 In the context of AC1S4 this transition is nothing more than the Hawking-Page transition [85]. 
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(5.4) 



CBHstab = 
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Implementing the constraint of putting the system in a box by hand is not too hard. For 



yo, and fixed total charge. We can then repeat the existence and stability analysis done for the 
asymptotically flat case. This is shown in fig. 5.2, the boxed analog of fig. 5.1 (not plotting the blue 
region). All the qualitative features remain intact, except that the phase boundaries are pushed 
down significantly. In fact, they are pushed down below the dotted line indicating the maximal 
temperature T as a function of c at fixed Qq. Thus, all bound states occur in the region of positive 
specific heat, and we can say that at sufficiently high temperature, there is only the black hole 
solution, and it is stable. For the infinite system we could only make the analogous statement for 
sufficiently high nonextremality parameter (mass) of the black hole. If we imagine the existence of 
some holographic dual field theory description of the system under consideration, with the black 
hole representing the unique disordered high temperature thermodynamic equilibrium state and 
the multiple bound state configurations as the (meta)stable thermodynamic states characteristic 
for glassy systems below their critical temperature, then a positive specific heat throughout the 
parameter regime of interest is certainly an expected feature. 

5.4 Thermodynamics and phase structure 

Figures 5.1 and 5.2 are certainly suggestive of an interesting phase structure, with the attractor 
points of the extremal black holes at (Qo,c) = (±g,0) corresponding to quantum critical points, 
the grey regions corresponding to the "normal" disordered high temperature phase and the colored 
regions to glass- like phases. 

However, to say anything definitive about phase structure of this system, a more careful study 
of the thermodynamic weight of various configurations is needed, as well as an analysis of candidate 
order parameters and how they scale near phase boundaries and near the critical points. We will 
leave a full analysis to future work and restrict ourselves here to some simple observations. 

At fixed temperature, a thermodynamic system will try to minimize its free energy F = E — TS. 
The free energy satisfies the first law 



where 5W denotes the work delivered to the system in some infinitesimal process during which the 
temperature changes by an infinitesimal amount ST. In particular this means that if we move a 
probe particle from some position r = R into the black hole, while keeping the temperature fixed, 
the total change in free energy of the system 5F = 5Fbu+SF p must be equal to 5W = Vp(horizon) — 




5F 



S5T + 5W 



(5.5) 
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V P (R), with V p = V g + V em the probe potential defined in (4.1). 7 Since V^(horizon) = 0, this means 

F(probe at R) - F(probe in BH) = V p (r = R) , (5.6) 

and thus at fixed temperature the configuration with the probe outside of the black hole will be 
thermodynamically favored over the one with the probe inside simply whenever V P {R) is negative. 

As a check, we note that for R = oo, this is the statement that 5Fbr + SF p = — V^(oo) — V era (oo). 
Since 5F P = —m p = — V g (oo) (as the probe is first at infinity and then gone, and we are allowed to 
ignore as before thermal contributions to the probe free energy), this is equivalent to SF^n = — V em , 
which can be directly checked from the expressions in 3.5, without using the first law, taking care 
to vary c at the same time as the charges to keep T fixed. 

To conclude, figures 5.1 and 5.2 give information about the thermodynamic preferences of the 
system for perturbations around the single black hole state. For example when in the green region, 
the black hole will start to populate black hole halos of many different charge types. We should keep 
in mind however that as soon as the number of such probes becomes macroscopic, or when they 
coalesce into black holes of size comparable to the background black hole, our neglect of the thermal 
internal, kinetic and inter-probe interaction energies is no longer justified. Hence we cannot read 
off the endpoint of this evolution from the diagrams. Nevertheless, the existence of exponentially 
many bound states with free energy below the black hole free energy strongly suggests glass-like 
behavior. This is not entirely obvious though, since the plots also show that the black hole core 
of any negative energy bound state is always unstable to emission of particles to infinity or to 
the boundary of the confining box, so for a sufficiently large box, it is conceivable that the final 
equilibrium state may still be a simple dilute charged gas. We will revisit such questions in future 
work. 

Finally, let us comment on the tentative interpretation of the points Qq = ±g at T = as 
quantum critical points. As we mentioned several times already, they correspond to the case in 
which the scalars at infinity are at the black hole attractor point, yo = y* = ^/ 6 |p^ . In this case, 
the background scalars do not flow but rather are constant over all of space, and no molecular 
bound state configurations exist - everything has been sucked into the black hole or pushed away 
to infinity. (It is known that this remains true away from the probe limit, at least in the BPS 
case: non-marginal bound states cannot exist at the attractor point of the total charge.) At zero 
temperature, the geometry is that of an extremal Reissner-Nordstrom black hole, developing in 
particular an AdS2 X S 2 throat, suggestive of a holographic conformal fixed point. 

Taking yo away from the attractor point but still keeping T = 0, the AdS2 is preserved but the 
scalars will now flow from yo a t infinity to at the horizon. Putting the system in a box of finite 
radius R and decreasing R while keeping yo fixed will have roughly the same effect as moving yo 

7 Recall we are treating the probe as a structureless object without internal thermal energy. See footnote 4 for the 
argument of why this is justified in the probe limit even if the probe itself is a black hole. The argument extends to 
thermal kinetic energy of the probe moving around in its potential well (provided this lies at finite rescaled radius), 
which we thus ignore as well. 
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along the attractor flow towards the fixed point y+; in particular near Qq = g, the configurational 
entropy will decrease by some power of R dictated by (5.3). In a hypothetical field theory dual, 
regardless of the large r asymptotics of spacetime [88], decreasing R would correspond to flowing 
to the IR, so this would correspond to a power law decrease of the configurational entropy when 
coarse graining over increasing length scales. 

If |<5o| = g is to be a quantum critical point associated to a quantum phase transition, then 
we should see different physical properties on the two sides of it. Near the supersymmetric critical 
point, one such property can be inferred by inspecting (5.1). Recall that pulling a probe out 
of the background black hole produces a D6 magnetic dipole moment pointing along the radial 
direction, equal to ji = po r eq = PO\/y~o~Pi -75 ^°~f 1 , as well as an angular momentum J = ^(7, T) = 

^PiviPoiQo - Hence 

J = gjl, g = ^y 3/2 (<?i - §) . (5.7) 

From (5.1) we see that the range of possible values of \g\ runs from to a maximal value proportional 
to I Qo — g I, and that it changes sign across the phase boundaries. The coefficient g is a physical 
observable, and could be used as an order parameter to distinguish the two putative phases. 
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A Counting configurations 

There exist supersymmetric probe bound state solutions for all values of Qo > except 1/6. 
The number of possible bound state solutions will not be constant however. In particular when 
Qq — > 1/6, the allowed region in the probe charge space shrinks to zero. In this appendix we obtain 
an estimate for the number of supersymmetric probe bound state solutions near this point. First 
we consider the case in which we simply count the number of allowed probe charges, later on we 
will include the lowest Landau level degeneracies for each choice of probe. 

A.l Single probe 

For a single probe, this number is A/i = X^eA 1> wnere A is the allowed region in probe charge 
space, bounded by the requirement that the bound state exists and the probe approximation is 



23 



satisfied. As we take the limit Pi — > oo, A will contain an increasingly large number of charges, 
hence the number of lattice points contained in A can be estimated by computing the volume of A 
in charge space: 



dp dpi dqi dq 1 



dpoPoVo I dbdhdk. 
A 



(A.l) 



Here A is the region of allowed values of {k,n, b). Since \k\ is bounded by y 2(Qo + b) and y 2(g + 6) 

according to (5.1), the integral over k gives a factor fi(b) ~ \/ 1 + b — \J Qo + b . (We drop 

irrelevant numerical factors here and in what follows.) Because of the constraint n 2 < |& 3 , the 
integral over h gives a factor /2(b) ~ 6 3//2 . Performing the integral over b gives 



rf6/i(&)/ 2 (6)~|Qo 



i| 5 2 

6 I "max 



(A.2) 



Here 6 max is the maximal value of b, which we have assumed to be very large so we are allowed 
to drop terms of order 6 max and lower. The reason 6 max is not infinite is the requirement that 



the probe approximation should be valid. To estimate it, recall that m p /M 



yopo ™>p 
Pi M 



Since we 



are exploring the region Qq 
approximation requires 6 max 



1/6, M is of order 1. For large b, fh p 



b 3 / 2 , so the probe 



'yoPo J ' ~' ' p 

maximal m p /M we allow. We are then left with the integral over pq: 



where a = 2/3 and e ~ m p /M is some small number, the 
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Pi 
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(A.3) 
(A.4) 
(A.5) 



The upper integral bound po < eP\ /yo comes again from requiring that we remain within the probe 
approximation, this time in the limit of large po. Notice that the final result does not depend on the 
actual value of a. In fact, we could have inferred the prefactor simply from the scaling symmetries 
of the system: From (A.l) it is clear that A/i has scaling weights (0,6,4) under the symmetries 
of 4.2, while e, Pi and yo have scaling weights (—1,0,1), (1,1,0) and (0,1,0), respectively. This 
uniquely determines their powers. 



The above expressions are valid when 

y* = 



Qc 



< 1, i.e. 



6Q0 
Pi 



On 
Pi ' 



close to the attractor point: yo — * V* = *\J^~pf~- Notice that since the probe approximation requires 
yo "C Pi (as discussed in section 4.2), self-consistency in this regime requires Qq <C Pf, i.e. we are 
necessarily in the non-Cardy regime. 

Pushing yo an d e as high up as possible while conceivably still yielding more or less sensible 
results, i.e. 

yo~Pi, Qo~Pi 3 , e~l, (A.6) 
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we get 

M.max ~ Pf ■ (A.7) 

A. 2 Multiple probes 

In the probe approximation we can also easily build multi-probe bound states: By assumption, the 
probes do not backreact so we can simply superimpose the single probe configurations, as long as 
we keep the total probe mass ^ m p i small compared to M. 8 Imagine a general situation in which 
the number of single probe states with m p /M < e is given by 

Af 1 (e) = Ae n , (A.8) 

where A is some large number. In the case at hand n = 4, but we will keep things slightly 
more general here for future reference. The density of single particle states at -»# = e is then 
dMi(e) = ne n ~ 1 A. For K labeled probes, the density of states at ^ = ej, i = 1, . . . , K is 

dAT K (ei, ...,e K ) = {nA)*^ 1 ■ ■ ■ e™" 1 de x ■ ■ ■ de K . (A.9) 

The total number of states with an arbitrary number K of unlabeled probes satisfying £\ ?jp < e 
is therefore 

The X/K\ corrects in a classical way for overcounting. 9 The integral can be factorized by represent- 
ing the constraint £^ e% < e as the contour integral ^ / j^e x ( e ~^ i6i \ where the contour is taken 
to be on the right of the pole at A = 0. This yields 



J K=0 yJ 

1 f d\ / n\A\ , k . 

exp Ae+— . (A.ll) 



2m J X V A 

At large A this can be computed by saddle point evaluation. To leading order, dropping order one 
numerical factors: 

logAA(e) ~ (Ae n )^ . (A. 12) 



8 In general the probes will interact with each other, with interaction strength given by their mutual symplectic 
products. For probes which happen to have small symplectic products with the background black hole, i.e. probes 
which are close to be swallowed by the black hole, these interactions become important even in the probe limit. We 
will ignore such boundary cases here. 

9 We do not use quantum statistics because the probability that two probes occupy the same quantum microstate 
is completely negligible in this setup. This is already true if the probes are considered to be point particles without 
internal degrees of freedom, but becomes obvious without work when one takes into account the huge number of 
internal microstates the individual probe black holes can choose from. 
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Applying this to (A. 5) gives for the total number of distinct probe configurations 

me)~exp( K e^\yl-yl\ 1/5 P? /5 ) , (A.13) 

where k is some order 1 constant and y* = is the attractor fixed point. 
In the regime (A. 6), we thus get a configurational entropy 

logAA max ~ Pf 75 . (A. 14) 

A. 3 Including Landau level degeneracies 

A single probe bound to the background black hole has classically an S 2 moduli space, but because 
a magnetic field threads the sphere, the space of quantum BPS ground states (i.e. the lowest Landau 
level) will be degenerate. The degeneracy is given by the effective magnetic flux as seen by the 
probe [13]: 

d 7 = \poQo - 01 -Pi | = PoPiVolQo ~qi\- (A. 15) 

To count the total number of such one particle ground states (ignoring the internal degrees of 
freedom of the probe and of the black hole), it suffices to replace the 1 in (A.l) by the LLL 
degeneracy factor eL. Because this has scaling dimensions (1, 3, 1), doing so will add an additional 
factor eyoP 2 to the final result (A. 5). Furthermore, because the insertion | Qo — <7i| is of generically 
of order |Qo — g| over the integration domain, its effect will be to modify the power of |Qo — g| 
from linear to quadratic. All in all we get 



2 



A/- LLLil (e) ~ e 5 Q - \ y 3 P? (A. 16) 

2 p6 

yo ' 



6Qo ,,2 



(A.17) 



For the multi-probe system we get from this, using (A. 12), 



log A/lll(c) ~ 6 5/6 ^ " yo ^ (A.18) 
Finally, in regime (A. 6), we get a configurational entropy 

logA/LLL,max ~ P\^ ■ (A. 19) 

Recall that the entropy of the black hole scales as Pf . 
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